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Abstract. We reanalyse the spectra used by D.-E. Lieb- 
scher et al. (1992) with the same goal - determining the 
cosmological parameters - and basically the same assump- 
tions but with a different statistical method, which does 
not rely on binning the data. Also, we correct for selection 
effects. We basically confirm their result, with somewhat 
larger but still very small (formal) errors. However, all 
world models within the 99% confidence region are ruled 
out because a firm lower limit on fig is significantly larger. 
This directly demonstrates for the first time the existence 
of intrinsic evolution in the Lyman-a forest. 
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1. Introduction 

A method of determining the cosmological parameters Ao 
and O has recently been suggested by D.-E. Liebscher 
et al. ( |1992j - hereafter LPH). They use the rcdshift distri- 
bution of Lyman-a forest lines, with the assumption that 
the Lyman-a forest absorbers have a constant comoving 
density. This basically determines the volume element as 
a function of redshift, which depends on the cosmological 
model and is especially sensitive to a positiv e valu e of the 
cosmological constant Ao- (See, e.g., Feige ( 1992 ) for an 
explanation of this effect.) Under this assumption, which 
means no intrinsic evolution - we henceforth use 'evolu- 
tion' to mean intrinsic physical evolution and not in the 
more general and somewhat misleading sense of 'depen- 
dence on z' - they find an acceptable fit only in a small 
area of the Ao-f2o~P iane - It is important to realise that, if 
evolution is present in the absorbers, one would not neces- 
sarily expect to find any (then necessarily wrong) cosmo- 
logical model (not even an extreme one perhaps excluded 



Send offprint requests to: P. Helbig 
* rkayser@hs.uni-hamburg.de 
** phelbig@hs.uni-hamburg.de 



by other arguments) which would be compatible with the 
observations without requiring evolution; i. e., evolution- 
ary effects and the influence of the cosmological model 
are not necessarily degenerate, especially when, as is here 
the case, one has a large data set extending to high red- 
shifts. Of course, the fact that LPH find an acceptable fit 
docs not prove the assumption of no evolution, but it does 
mean that one cannot say that 'evidence of evolution is 
directly seen in the spectra'; the fact that an acceptable 
fit is found demonstrates that a cosmological model exists 
in which no evolution is required to explain the observed 
spectra. 

The paper by LPH was subjected to much criticism, 
but our investigation shows that most of this is unfounded. 
We have confirmed their results with a better statisti- 
cal method (and have reproduced their results with their 
method as a check at the outset of our investigation) ; the 
result is perfectly consistent with the assumptions and 
gives an acceptable fit. 

This was admittedly somewhat surprising, since we 
have a number of points of criticism of LPH. First, the dis- 
tribution of spectral lines was not used directly, but rather 
the lines were counted in redshift bins. Binning introduces 
free parameters and throws away information (cf., e.g., 
Press et al. ( |1992| )) and should be avoided if possible. Also 
on the mathematical side, a fit to their Z 2 instead of Z 
was done. In addition, since the smallest equivalent width 
capable of being registered by all detectors used in the 
inhomogeneous sample of Lyman-a forest lines is 0.2 A, 
this implies, because the lowest redshift in their sample 
is z = 0.003, that the intrinsic equivalent width H^ntr to 
which the sample is complete must also be w 0.2 A. Of the 
1297 lines used by LPH only 832 have a W- m tr larger than 
this minimum. Furthermore, since the physical distance 
between two absorbers is a function of their redshifts and 
the cosmological parameters Ao and f2o, their sample is 
only complete to distances greater than that correspond- 
ing to the minimum observable distance in the given world 
model. To avoid biasing, we must for each world model ex- 
clude a certain number of lines from the statistical anal- 
ysis. (See the di scuss ion of the correction for 'resolution 
bias' in Kayser fll995| ).) 
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2. Basic theory 

2.1. Cosmology 

We make the 'standard assumption' that the Universe can 
be described by the Robertson- Walker metric. With the 
assumption of no evolution, the absorbers have a con- 
stant comoving density. The prope r distance (this is the 
distance D p in Kayser et al. ( 1996 ), hereafter KHS) in a 
given world model between two absorbers with redshifts 
z x and z y is given by D p (z y ) — D p (z x ). This is simply 



D p = 

xy 



where 



Q(z) = fi„(l + zf - (Oq + Ao - 1)(1 + zf + A . 



(1) 



(2) 



The idea is to find the cosmological parameters for which 
D P y corresponds most closely to what is observed in the 
spectra; if there is no evolution, then this must be the 
correct world model, if the fit is acceptable. We use the 
distance between adjacent absorbers Az as the observa- 
tional quantity. This is of course equivalent to the number 
of lines per redshift interval except for the loss of precision 
due to binning in the latter case. 



2.2. Statistics 

We use a method similar to that described 



Kayser (1995): We calculate the distribution of distances 
between pairs of adjacent absorbers, divide the sam- 
ple into N subsamples based on redshift (with the 
same number of objects in each subsample) and com- 
pare m — N/2 independent pairs of samples with m 
Kolmogorov-Smirnov (K-S) tests (instead of just two sam- 
ples and one test as in Kayser ( |l995| )n) . Since the null hy- 
pothesis is that the distribution is independent of redshift, 
the choice of independent pairs is arbitrary. 

What is the probability P m i n that the minimum of all 
m K-S probabilities Pj is less than a given Pq? Since, if 
min (P^ > Pq all values Pi must be larger than Pq we 
have 



P(min(P) > P ) = (1 - P ) m 
and thus 

P min = P(min(P) < P ) = 1 - (1 - PoY 
Solving for Pq we obtain 

PQ = 1 - Vl " P mi„ 



(3) 



(4) 



(5) 



1 It is obvious that using just two subsamples will not allow 
one to reject the null hypothesis that the distribution is inde- 
pendent of redshift if there is a symmetry with respect to the 
redshift separating the subsamples, even though the distribu- 
tion is not independent of redshift. 
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Fig. 1. Demonstration of the statistical method. Plotted is 
the comoving distance Y in units of c/Ho between pairs of 
adjacent absorbers as a function of redshift. Also shown is the 
minimum observable distance (lower curve) in the given world 
model in the interval where data points exist; the maximum of 
this curve (in the interval where data are used) determines the 
(straight horizontal) line below which the data cannot be used 
if one requires a complete sample. Here, Ao = and Qo — 1. 
Pairs of subsamples in different redshift intervals are compared 
with Kolmogorov-Smirnov tests and a probability is calculated 
that the different subsamples could have been drawn from the 
same parent population, thus giving the probability that the 
distribution is independent of redshift 



We can reject the null hypotheses that this distribution is 
independent of redshift at the 100(1 — P m in)% confidence 
level in the cases where min(P^) < Pq. Here, we conserva- 
tively considered a probability P min of 1% or larger to be 
compatible with the observations and used three pairs of 
samples, thus ruling out all world models with a value of 
min(P) less than P = 1 - v^I99 » 0.0033 at the 99% 
confidence level. Of course, as mentioned above, for each 
world model the actual sample is different in order that it 
be complete down to a certain distance between pairs of 
adjacent absorbers. (See Fig. |l|.) 

3. The sample 

For purposes of comparison, we used the same sample of 
spectral lines as in LPH, except that we omitted 7 lines 
with z < 1.65 (taken from Bahcall et al. ( |l99l| )). This 
means that the sample is now complete to 



Wi n t r = 0.2/(1 + 1.65) « 0.0755^ 
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Table 1. Spectra used. The redshift z refers to the midp oint o f 
the wavelengt h int erval. References are Atwood e t al. ( 1985), 



Pettini et al . ([19901 ), Sargent et al. ( [1988|), Steidel ( J1990[ ), Car- 
swell et al. Ql99ip and Rauch et al. (|1992j) 



# 


QSO 


z 


A [A] 


AA 

A 


reference 


1 


0420-388 


2.25 


3850- 


4050 


0.5 


Atwood 


2 


0420-388 


2.46 


4100- 


4300 


0.5 


Atwood 


3 


0420-388 


2.46 


4450- 


4650 


0.5 


Atwood 


4 


0420-388 


2.99 


4750- 


4950 


0.5 


Atwood 


5 


2206- 199N 


2.25 


3860- 


4052 


0.09 


Pettini 


6 


2206- 199N 


2.43 


4070- 


4282 


0.09 


Pettini 


7 


0114-089 


3.11 


4900- 


5100 


1.5 


Sargent 


8 


0913+072 


2.45 


3900- 


4500 


1.5 


Sargent 


9 


1159+124 


3.11 


4600- 


5400 


1.5 


Sargent 


10 


1247+267 


1.84 


3260- 


3660 


0.8 


Sargent 


11 


1511+091 


2.70 


4300- 


4700 


1.5 


Sargent 


12 


1623+269 


2.29 


3700- 


4300 


0.8 


Sargent 


13 


2126-158 


3.11 


4800- 


5200 


1.5 


Sargent 


14 


0142-100 


2.45 


3900- 


4500 


0.8 


Sargent 


15 


0237-233 


2.04 


3500- 


3900 


0.8 


Sargent 


16 


0424-131 


2.04 


3600- 


3800 


0.8 


Sargent 


17 


1017+280 


1.80 


3300- 


3500 


0.8 


Sargent 


18 


2000-330 


3.44 


5100- 


5700 


1.35 


Steidel 


19 


0055-209 


3.52 


5400- 


5600 


1.35 


Steidel 


20 


0000-263 


3.77 


5400- 


6200 


1.35 


Steidel 


21 


1208+101 


3.60 


5400- 


5800 


1.51 


Steidel 


22 


1100-264 


1.91 


3440- 


3640 


0.11 


Carswell 


23 


1100-264 


2.05 


3640- 


3780 


0.11 


Carswell 


24 


0014+813 


2.74 


4500- 


4600 


0.36 


Rauch 


25 


0014+813 


2.87 


4600- 


4800 


0.36 


Rauch 


26 


0014+813 


3.03 


4800- 


5000 


0.36 


Rauch 


27 


0014+813 


3.15 


5000- 


5100 


0.36 


Rauch 


28 


0014+813 


3.28 


5100- 


5300 


0.36 


Rauch 



which gives us a sample consisting of 1108 of the 1297 LPH 
lines; including these 7 low-redshift lines would decrease 
the complete sample, as explained in Sect. ||, to 832 lines. 
Although the low-redshift lines have a significant impact 
on the x 2 method used by LPH, leaving out this few lines 
makes practically no difference in the statistic we use, so 
we choose to leave them out in order to be able to work 
with a larger complete sample. Since we wanted to use the 
actual redshifts of the lines and not simply the number 
of lines per redshift interval, the wavelengths of the lines 
were determined from scans of the literature spectra and 
the redshifts (and Wi ntr ) were then calculated. See Table [l] 
for more details. Our measured wavelengths of all the lines 
are available from the following URLs: 



http:/ /www. hs.uni-hamburg.de/cnglish/pcrsons/hclbig/ 
Research /Publications/Into/ ly man-alpha, html 



ftp:// ftp.uni-hamburg.de/pub/misc/astronomy/ 
lyman-alpha.tar.gz 
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Fig. 2. The probability as a function of the cosmological pa- 
rameters. The (thick) contour line is at 1%. Only world mod- 
els within this area are compatible with the data. Note the 
extremely small area of the Ao-fio-plane in this plot - out- 
side this area, no areas with an acceptable probability exist. 
The lower left corner of the plot corresponds to the de Sitter 
model. 2 For comparison, the LPH 1 a and 3 a error ellipsoids 
are included 



4. Calculations, results and discussion 

First, for a general overview, we calculated the probability 
for cosmological models in the area given by — 10 < Ao < 7 
and < tt < 7, with a resolution of 0.06 (« 33 000 dif- 
ferent models). Almost the entire area of this plane is 
rejected. Then, in a series of steps, the area containing 
a probability P m in > 1% was examined in progressively 
higher resolution. Only a small region, shown in Fig. ^ 
was found to be compatible with the observational data. 

Thus, we find in agreement with LPH that - with the 
assumption of no evolution cosmological models exist 
which provide an acceptable fit to the data. The fact that 
the area in the Ao-^o-plaue in which this is the case is so 
small means that, if the method is correct, it offers by far 



2 Most of the space to the lower right of the allowed area 
corresponds to world models in which the maximum possible 
redshift is less than the largest Lyman-a redshift in the sample. 
Most of these are the so-called bounce models, which have no 
big bang and, for realistic values of Qo, are excluded by the 
mere existence of obje cts w ith a large cosmological redshift 
(see Stabell & Refsdal ( |l966| ) or KHS). For low fi values there 
are some bounce models at the lower part of the allowed area 
which are not excluded by the redshift of the lines used in this 
analysis. 
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the best known method of determining the cosmological 
parameters. The fact that the allowed area is somewhat 
outside the (present) canonical range of cosmological mod- 
els is probably one reason for some of the criticism of the 
work of LPH. Mostly, this has concentrated on the fact 
that this cosmological model is dominated by the cosmo- 
logical constant. We hasten to point out that it has not 
yet been shown that there is a way of excluding this large 
a value of A in spatially closed (k = +1) world modelsE 
It is generally recognised that gravitational lensing statis- 
tics can provide the tightest constraints on a positive cos- 
mological constant. However, at present such constraints 
have only been d erived for flat (k = 0) models. (See, e.g., 
Kochanek ( 1996 ) and references therein.) Basically, large 
Ao values predict too many lenses, primarily because of 
the relatively large volume element dV/dz in these mod- 
els. For k = +1 models, this effect is not as pronounced, 
since dV/dz is generally smaller in a k = +1 model than 
in a flat model with the same value of Ao. Thus, a given 
upper limit on Ao for a flat model corresponds to a weaker 
limit for a k — +1 model. Nevertheless, it is.-conceivable 
that lensing statistics could rule out the LPHcl model and 
other k = +1 models with a large cosmological constant; 
this is currently under investigation. 

LPH point out 'advantages' of their model, such as the 
fact that a relatively old universe and a relatively large 
Hubble constant can both be accommodated. In fact, bet- 
ter values for the world age and H will probably rule out 
some standard cosmological models in the near future. The 
point is that such indirect tests don't rule out the LPH 
model; if anything, they point in the direction of this or a 
similar cosmological model. 

There is not necessarily any non-baryonic matter in the 
LPH model, since their value of Qq, coupled with plausible 
values for Hq, gives a density comparable to the baryonic 
densit y predicted by nucleosynthesisH (See, e.g., Krauss 
( 1995 ) and references therein.) This could also be seen as 
a point in favour of this cosmological model. 

However, the « 3tr upper limit for the total density Qq 
(« 0.035) is significantly less than the extremely conser- 
vative lower limit (w 0.05) advocated by Peebles ( |1993| , 
Chapter 20). Since this conservative lower limit on f^o, 
based on well-understood physics, is considerably less than 
lower limits advocated by other researchers, this means 
that the LPH model is definitively ruled out and that it 



3 All cosmological models in Fig. ^| have k = +1 and Ao > 1 
and thus are spatially closed and will expand forever. See KHS. 

4 In the following, the term 'LPH model' should be under- 
stood to mean the range of cosmological models within the 
errors of the LPH best fit model as well as our results since, 
for the purposes of the following discussion, the cosmological 
models are identical. 

5 For some 'more standard' low-density standard cosmolog- 
ical models, it is also just barely possible to have no non- 
baryonic matter, but only for an extremely low value of Ho 
and a baryonic density at the upper limit. 




Fig. 3. Evolution in the Lyman-a forest for the Ein- 
stein-de Sitter model (Ao = 0, Qo = 1). Higher values mean 
greater distances between absorbers and thus fewer absorbers 
per comoving volume. (Values are normalised with respect to 
the mean.) Different curves correspond to different binning 
schemes 



is not possible to determine the cosmological parameters 
by this method. 

5. Evolution 

The conclusion in the last paragraph means that there 
must be evolution in the Lyman-cv forest, because without 
evolution, the only cosmological model compatible with 
the observations is definitively ruled out on other grounds. 
This conclusion does not assume any a priori values for 
the cosmological parameters. 

For a given cosmological model, the comoving dis- 
tances between adjacent Lyman-cv forest absorbers can be 
calculated, as shown above. A plot of the average distance 
per redshift bin should be almost flat if the cosmologi- 
cal model used for the calculation is correct and if there 
is no evolution. Previously, we assumed no evolution to 
so determine the cosmological model. However, since we 
now know that there must be evolution in the Lyman- 
cv forest, we now, for a few representative cosmological 
models, calculate what evolution is necessary to explain 
the observations. Several different binning schemes were 
used to isolate real effects from noise. (See Fig. |). The 
resulting calculated evolution is relatively independent of 
the world model (except near the LPH model, where of 
course the calculated evolution is almost flat). This means 
that Fig. |^ is approximately correct for almost any world 
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the A 2 curveO (following Stabell & Refsdal fll966| )) are the 
so-called plateau models, where the universe has a quasi- 



Fig. 4. The same as Fig. |§] but with one binning scheme and 
different curves correspond to different cosmological models. 
The thick curves correspond to (Ao,f2o) values of (0.0, f.O) 
(0.0,0.3), (0.65,0.35) and (2.0,4.0). The curves are practi- 
cally the same within the binning errors (cf. Fig ^). The thin 
curve corresponds to the best-fit model of HPL, (Ao,f2o) = 
(f. 080,0. Of 4), for which of course the calculated evolution curve 
is almost flat within the binning errors. 



model (cf. Fig. [|). This is the first calculation of the evo- 
lution of the Lyman-a forest based on observational data 
which makes no assumptions (other than basic cosmologi- 
cal ones). The general trend is an increase in the comoving 
density of Lyman-a forest absorbers with z and thus a de- 
crease with time. 

The form of the curves in Fig. || make it possible to 
understand why the LPH model is compatible with the 
observations without requiring any evolution. The basic 
feature in the curves in Fig. || is the fact that the comov- 
ing distance between adjacent absorbers decreases with 
increasing redshift. If there were no evolution, then the co- 
moving distance should be independent of redshift. Thus, 
a decrease in the comoving distance between adjacent ab- 
sorbers with increasing redshift can be compensated in a 
world model in which the comoving distance for a fixed 
Az is comparatively larger at larger redshifts. This is only 
the case in a small part of the Ao-f2o-pl ane - 

For a given Q Q value there is a value of Ao above which 
no big bang is possible; rather, the universe contracts 
from oo to a finite R m in before expanding again. Near 
this curve separating world models with and without a 
big bang, which for brevity we refer to in the following as 



static phase. (See, e.g. Feige (1992) for a discussion of this 
phase.) Near the redshift of the plateau, the comoving dis- 
tance for a fixed Az can become very large. The redshift 
of the plateau depends sensitively on the value of the cos- 
mological parameters, ranging from oo near the de Sitter 
model (Ao = 1, = 0) to values < 1 for relatively small 
values of U and Ao- Figure || shows the comoving dis- 
tances between adjacent absorbers becoming smaller with 
increasing redshift, here assumed to be due to intrinsic 
evolution, ft is easy to see that a cosmological model - 
such as the ones near the LPH model - in which the dis- 
tance is relatively larger at larger redshift will counteract 
this effect, producing a flatter distribution which, in such 
a cosmological model, would imply a lack of intrinsic evo- 
lution. Since the redshift of the plateau and hence of the 
large volume element can be adjusted by moving along the 
Ai curve, and its width by moving perpendicular to it, it 
is easy to find a cosmological model which, with no evo- 
lution, will be compatible with observational data from a 
'normal' cosmological model in which the comoving dis- 
tance between absorbers for a fixed Az decreases with 
increasing redshift due to evolutionary effects. Due to the 
sensitivity mentioned above, this will only be possible in 
a very small area of the parameter space. Since the width 
of the plateau is even more sensitive to the distance from 
the Ai curve than is its redshift to the position along the 
curve, the allowed area in Fig. || has an elongated shape. 

6. Summary and conclusions 



We confirm the result of D.-E. Liebscher et al. (1992): 
with the assumption of no evolution, observations of the 
Lyman-a forest allow only a small range of (spatially 
closed) cosmological models with a positive cosmological 
constant and a low value of £}q. This is despite the fact 
that we avoided some of the pitfalls of their method such 
as binning and selection effects. The former is relatively 
harmless in this case because of the large number of data 
points; the fact that the latter isn't important is probably 
just coincidence. 

However, the fact that all cosmological models within 
the errors of the 'best fit' model are ruled out because a 
firm lower limit on flo is significantly larger means that 
there must be intrinsic evolution in the Lyman-a forest. 
Calculating the necessary evolution for a variety of cos- 
mological models gives similar results, so that one can 
make some qualitative conclusions about the Lyman-a 
forest with almost no knowledge about the cosmological 
model: the comoving distance between adjacent absorbers 
decreases with redshift (and thus increases with time). 
This is a result which is very robust in that, in the area of 



6 ThisisWMTYPE 17 as given by the INIC0S routine mentioned 
in KHS. 
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parameter space which is not excluded on other grounds, 
it is almost independent of the cosmological model. 
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